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Abstract 

This work examines a novel heterogeneous dynamic spectrum access network where the primary 
users (PUs) coexist with both underlay and interweave cognitive radios (ICRs); all terminals being 
q . potentially equipped with multiple antennas. Underlay cognitive transmitters (UCTs) are allowed to 

transmit concurrently with PUs subject to interference constraints, while the ICRs employ spectrum 
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£> ■ sensing and are permitted to access the shared spectrum only when both PUs and UCTs are absent 
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p — I . cognitive receivers (UCRs) and constraining interference leaked to PUs. The objective of such a proactive 

approach, referred to as prescient precoding, is to minimize the probability of interference from ICRs 



o 



X 



We investigate the design of MIMO precoding algorithms for the UCT that increase the detection 
probability at the ICRs, while simultaneously meeting a desired Quality-of-Service target to the underlay 



to the UCRs and primary receivers due to imperfect spectrum sensing. We begin with downlink 
prescient precoding algorithms for multiple single-antenna UCRs and multi-antenna PUs/ICRs. We 
then present prescient block-diagonalization algorithms for the MIMO underlay downlink where spatial 
multiplexing is performed for a plurality of multi-antenna UCRs. Numerical experiments demonstrate 
that prescient precoding by UCTs provides a pronounced performance gain compared to conventional 
underlay precoding strategies. 
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I. INTRODUCTION 

Dynamic spectrum access (DSA) is emerging as a promising solution to enable better utiliza- 
tion of the radio spectrum, especially in bands that are currently under-utilized Q. DSA partitions 
wireless terminals into categories of primary (licensed) and secondary (cognitive radio) users, 
where the primary users have priority in accessing the shared spectrum. Furthermore, the two 
most prevalent classifications of secondary users are underlay cognitive radios and interweave 
cognitive radios (ICRs), following the terminology of 0. The underlay paradigm mandates that 
concurrent secondary and primary transmissions may occur only if the interference generated 
by the underlay cognitive transmitters (UCTs) at the primary receivers (PRs) is below some 
acceptable threshold. In contrast, ICRs are allowed to opportunistically use the spectrum only 
when it is not occupied by a primary transmitter (PT) with priority. In the absence of standard 
control channels or coordinated medium access between the primary and secondary users, the 
ICRs must periodically sense the spectrum for the presence of PTs [DOl-GDl and cease transmission 
upon detection. Inevitably, imperfect ICR spectrum sensing due to channel fading and other 
impairments will lead to unintentional interference at the underlay cognitive receivers (UCRs) 
and PRs. 

Underlay and ICR networks have been studied separately in extensive detail for both single- 
antenna and multiple-input-multiple-output (MIMO) terminals HI, 0. The use of multiple 
antennas in ICRs has been suggested for improved spectrum sensing capabilities by means of 
receive diversity Hl- lfTOll . MEMO systems have also been investigated in the context of underlay 
DSA networks, where multiple transmit antennas are used by UCTs for beamforming and to 
control the interference to the PRs, assuming either complete or partial channel state information 
(CSI) at the SU transmitter flTTJ - [JT5]] . However, there is little if any prior work that examines 
heterogeneous DSA networks with both UCTs and ICRs attempting to coexist simultaneously 
with primary users. Note that such a scenario is significantly different from hybrid secondary 
users that are capable of both underlay and interweave cognition lfT6l . ifPTll . 

Therefore, this work examines a fundamentally novel heterogeneous DSA network where the 
primary users share their spectrum with both UCRs and ICRs; all terminals being potentially 
equipped with multiple antennas. Specifically, we investigate the design of MIMO precod- 
ing algorithms for a underlay downlink network with multiple UCRs and interweave radios. 



The heterogeneous DSA network presents a myriad of conflicting objectives for the underlay 
transmitter, since it must mitigate the multi-user interference among its own UCRs, constrain 
the interference leaked to PRs, and ensure that the detection probability of the ICRs is high 
so as to preemptively avoid interference from them. Consequently, this paper is devoted to 
the design of novel precoding algorithms, collectively referred to as prescient precoding, that 
balance these competing objectives. The aim of prescient precoding is to reduce the probability of 
interference due to imperfect spectrum sensing from ICRs to the underlay and primary receivers, 
while simultaneously meeting their QoS/interference temperature requirements. Therefore, our 
contributions include the following: 

• We introduce a novel heterogeneous DSA network with both underlay and interweave radios, 
and demonstrate that existing underlay precoding schemes are highly suboptimal for such 
networks. 

• We design new prescient downlink precoding schemes of varying complexity for multiple 
single-antenna UCRs and MIMO PRs/ICRs. 

• We present prescient block-diagonalization algorithms for the MIMO UCR downlink chan- 
nel when all UCRs are equipped with multiple antennas in addition to the PRs and ICRs. 

. We demonstrate via numerical simulations that diverting resources from the UCRs to 
increase the detection probability at the ICRs can significantly suppress unintentional ICR 
interference. 

The paper is organized as follows. Section HI] introduces the mathematical model of the 
DSA network and the spectrum sensing performance of the ICRs. Prescient downlink precoding 
algorithms for the case of single-antenna underlay receivers are proposed in Section Unl Sec- 
tion |IV] outlines a prescient block-diagonalization algorithm for a MIMO downlink channel with 
multi-antenna underlay receivers. Selected numerical examples are shown in Section |Vj and we 
conclude in Section [VH 

Notation: We will use CA/"(0, Z) to denote a circularly symmetric complex Gaussian distri- 
bution with zero mean and covariance matrix Z, £{■} to denote expectation, vec(-) the matrix 
column stacking operator, (-) T the transpose, (-) H the Hermitian transpose, (■)^ 1 the matrix 
inverse, Tr(-) the trace operator, |-| the matrix determinant, and I an identity matrix of appropriate 
dimension. 



II. Mathematical Model 
A. Signal and Network Model 

Consider a downlink network with a t u -antenna UCT, K u single- antenna UCRs as its intended 
destinations where K u < t u , K multi-antenna ICRs with 77 antennas each, and a single PT-PR 
pair with t p and r p antennas respectively, as depicted in Fig. [T] The scenario with multi-antenna 
UCRs is presented in Sec. [IV] Multiple PRs can be accommodated in the current model by 
aggregating them into a single virtual PR. 

Assuming linear precoding, the UCT downlink transmit signal at time index t is written as 

I<u 

x(t)=^w A ,,(i)=Ws M (t), (1) 

fe=i 

where W e 0" xX " = wi . . . wk u is the precoding matrix whose columns represent 
individual beamforming vectors, and s u (t) E C A " xl is the collection of i.i.d. underlay in- 
formation symbols drawn from an M-ary phase-shift keying constellation with second-order 
statistics E {s u s^} = I. A power constraint Tr (WW H ) < P is imposed on the UCT signal. 
Furthermore, the UCT designs its transmit signal so as to ensure that the detection probability at 
the ICRs is satisfactory and the interference temperature at the PR remains below a pre-specified 
threshold £ p , as explained in Sec. flUl 

Suppressing the time index, the received signal at UCR k in the absence of ICR interference 
(i.e., with perfect spectrum sensing) is 

y k = h fc Ws u + n k (2) 

where h^ G C lxt " is the corresponding complex channel vector from the UCT, and n k is a 
circularly symmetric zero-mean complex Gaussian noise sample with variance a\ which includes 
interference from the PT. 

We list below the major assumptions regarding the heterogeneous DSA network. 
• We assume a partial CSI model at the UCT, which is defined to mean that the UCT always 
has knowledge of the instantaneous realizations of all the downlink channels ({h*.}^) and 
UCT-PR ({hfc} fc=1 ) channels, but may know only the distribution of its channels to the ICRs 
and the ICR-to-UCR channels. 



• The UCT has knowledge of the ICR transmit powers and the parameters of the spectrum 
sensing scheme deployed at the ICRs, which in practice are likely to be pre-defined by 
spectrum regulatory agencies. 

• There is no coordination between the UCT and ICRs. The UCT and PRs have a limited 
coordination in order to exchange CSI and tolerable interference limits. 

• All ICRs are half-duplex, which precludes for example simultaneous data transmission and 
spectrum sensing. We only consider in-band spectrum sensing; i.e., sensing is conducted 
on the same band that is used for data transmission. 

• The UCRs employ single-user decoding and treat all ICR/PT interference as noise. The 
interference from the ICRs is assumed to be instantaneous, i.e., the processing delay due 
to spectrum sensing is neglected. 

B. ICR Spectrum Sensing 

We assume that prior to transmitting, each ICR uses a binary hypothesis test based on Mri 
spatio-temporal samples to determine whether or not the band is occupied. Under the two 
hypotheses, the signal received by the i th ICR at time n is: 

H : Zi [n] = m.i [n] , n = 0, . . . , M — 1 (3a) 

Hi : Zi[n] — C[i[n] + m.i[n] , n — 0, . . . , M — 1 (3b) 

where irij[n] ~ CJ\f(0, efl) is temporally uncorrected background Gaussian noise of known 
variance and qi[n] represents the presence of a signal in the band. The Mrj complex samples 
are composed of 2Mrj independent real and imaginary components [20J. We are only interested 
in the case where qj[n] = {^[n], . . . , ^^[n]} is due to the UCT and also possibly the PT, in 
which case 

qij[n] = fi,jWs u [n] + d^-Sp [n] , j = 1, . . . , 77, (4) 

where channels fj^- E C lxt " from the UCT and djj G C lx4p from the PT are assumed to be 
invariant over the M samples, and s p G C ipXl is the PT transmit signal with total power P t . 

A broad range of spectrum sensing algorithms with varying levels of complexity and requisite 
a priori information have been proposed in the literature ll4l- lfT0l . On one hand, the optimal 
matched-filter detector has the most prohibitive requirements for CSI and PT signal information, 



while on the other hand non-coherent energy detection is the simplest possible detector since it 
only requires an accurate estimate of the noise variance ef. A range of composite generalized 
likelihood ratio tests (GLRT) and feature detectors lie in between these extremes. For this work, 
we assume the ICRs employ non-coherent energy detection due to its simplicity and the fact 
that it is unnecessary to distinguish between the UCT and PT signals. 
The test statistic and threshold test for the energy detector is given by JH 

M-i n n 

^ = E E \ Z{ j t n ]i 2; Ti ^ ^ ^ 

n=0 i=l Ha 

where A« is the detection threshold. We begin our development by analyzing the detection prob- 
ability Po,i at ICR i assuming deterministic channels and signals from the UCT and PT. Under 
the null hypothesis Ho, we see from (l3al) that z^j [n] ~ CM (0, ef), whereas under the alternative 
hypothesis H\ we have z it j [n] ~ CM (/J>i,j [n] = fjjWs u [n] + d it jS p [n], ef). Therefore, the test 
statistic Ti is the sum of the squares of 2Mrj independent real Gaussian random variables and 
thus follows a %-squared distribution under both hypotheses: 

Ti ~ix 2 2M ri undergo (6) 

T * ~ e h' 2 2 M ri (p) under ^i 

where the noncentrality parameter p = €^ 2 Y^ n =o Y7/=i \lM,j [ n ]\ is a function of W. 

Since we have an even number of samples 2Mrj (real and imaginary components of each 
sample), the false alarm probability follows immediately from the central chi-square cdf as |[22| 

p ^ = e ? E h (?) ' (7) 



where Aj = e^Q I (PA is chosen to satisfy a target false alarm rate Pf, and Q 2 (•) is the 

X 2Mrj ' X 2Mrj 

complementary cdf of the central chi-square distribution. The detection probability is given by 




PD,i = QMrrlVP;\Hr), (8) 



where Qm Ti (-, ■) is the generalized Marcum Q-function Il25i As the number of samples M 
grows large, T approaches a Gaussian random variable in distribution by the central limit 
theorem (CLT). Under hypothesis Hi the CLT yields T t ~ M (e?Mn + p, efMn + 2efp\ and 



the corresponding detection probability 

P Dti ~ Q l P | , (9) 

\e i y/e 2 Mr I + 2p t 

where Q (•) is the Gaussian Q-function. 

C. ICR Performance Prediction at UCT 

The ability of the UCT to predict the spectrum- sensing performance of the ICRs is an important 
ingredient of the prescient precoding paradigm. Under the partial CSI assumption, it is highly 
unlikely that the UCT has knowledge of the PT-to-ICR channel realizations and signals required 
to compute © or ©. A more plausible scenario is that the UCT knows the realizations of its 
channels {Fj} to the ICRs, and assumes the PT-to-ICR channels undergo Rayleigh fading with 
distribution d^ ~ CM (0, a 2 di l) Vi, j. 

Going one step further, the UCT may not have knowledge of the instantaneous realizations of 
its channels to the ICRs either. In order to gauge the energy detection performance of the ICRs, 
the UCT assumes a Rayleigh fading scenario such that f 4j - ~ CM (0, cr^l) Vi,j, and djj ~ 
CM (0,a di i) as before. Furthermore, the UCT and PT signals are each assumed to be drawn 
with uniform probability from a complex M-ary PSK constellation, and all channels, signals, and 
AWGN samples are mutually independent. Given these assumptions, the ICR samples Zij [n] arc 
distributed as independent Gaussian random variables [51 for both hypotheses. The false-alarm 
rate is clearly identical to that in © since it is channel-independent. Under %x, E {zij [n]} =0 
and a 2 z>i = var{^[n]} = 2<rJ ji Tr (WW H ) + 2P t t p a 2 - + e 2 . Thus, T,~ ^xifrr, and the 
corresponding average detection probability is 

^-^E^?' <10) 

r=0 \ z,i / 

From the UCT's perspective, a missed detection (Type II error) at any of the ICRs leads to 
interference at the underlay receivers, and this phenomenon plays a pivotal role in the prescient 
precoding principle. It will be useful to define the Bernoulli-distributed indicator function F t as 

1 with probability (1 - P D)i ) 
with probability (Pd,j). 



Fi therefore models the likelihood that ICR i unintentionally causes interference to the underlay 
and primary receivers, and is a function of W via Pd,%- 

Clearly, it is in the UCT's best interest to ensure that the probability of missed detection at 
the ICRs is made as small as possible, or equivalently, that the probability of detection is made 
as large as possible. To this end, we introduce the paradigm of prescient precoding in the next 
section in order to improve the reliability of the underlay downlink. 

III. Prescient Downlink Precoding 

It has been elegantly established that the capacity region of a conventional non-cognitive 
multi-antenna downlink channel without structured interference is achieved through non-linear 
dirty -paper coding, since all transmitted signals are known non-causally to the transmitter [|29l . 
However, linear precoding schemes for the multiuser downlink have been extensively studied 
due to their significantly lower complexity and near-capacity performance in certain regimes, 
and thus we focus on linear transmit preprocessing at the UCT. By definition, the UCT must 
limit the (instantaneous or average) interference it causes to the PR to a predefined threshold £ p : 

Tr (NWW ff N ff ) < £ p (12) 

if the instantaneous channel N G C TpXtu to the PR is known. 

The signal at an arbitrary UCR inclusive of ICR interference due to missed detections can be 
written as 

K u K 

yk = h k w k s Utk + } y hkWjSuj +y^i^Vfc,iSj ;i + n fc , k = 1, . . .,K U , (13) 

j^k i=l 

V v ' * 

intra - UCR interference ICR interference 

where v*^ ~ CM (0, cr^J.) and s/j G C r/Xl represent the (1 x 77) interfering channel and signal 
vector of power Pi from ICR i. We are interested in the characteristics of the aggregate ICR 
interference power at the k th UCR, defined as 

K 

4(W) = ^F i ||v fc ,|| 2 P l . (14) 

Taking the expectation of the ICR interference power in (fT4)) with respect to indicator functions 
{Fi}f =1 and the ICR-UCR channels {v k ,i\f=i y ields 

K 

h (W) = Y, (1 " Po,i) Pin°li- (15) 

i=l 



The UCR SINR that can be computed at the UCT is then approximated as 

IhfcwJ 
E^fclhfcW.I +4(W) + a| 

where the aggregate ICR interference /& (W) is a function of W via the spectrum-sensing 
detection probabilities. 

In the remainder of this section, we present several prescient design solutions for W that 
provide a tradeoff between complexity and underlay downlink performance. The attribute of 
"prescience" derives from the fact that the UCT anticipates interference at the PR from SUs due 
to imperfect spectrum sensing and takes preemptive measures to avoid the same. 

A. Direct UCR Sum Rate Maximization 

A wide variety of choices for W for conventional non-cognitive and underlay-only downlink 
channels have been explored in the literature. For example, a naive transmission scheme that 
disregards ICR CSI and PR interference would be to apply a modified regularized channel 
inversion (RCI) precoder 11301 , with 

Wo; = ^K (HuH? + K±l\ (17) 



where H„ = 



T 
n l • • • n K u 



v^ u V u P 

T 



, given the scale factor ( which is chosen as the smaller of the 
two scaling factors required to preserve the UCT transmit power and PR interference temperature 
constraints. However, the naive RCI precoder does not account for the potential ICR interference 
Ik, which can severely degrade the underlay sum-rate performance when I k is the dominant term 
of the denominator in (TT6l) . A more efficient usage of the side information available to the UCT 
is a direct sum-rate maximization approach that exploits knowledge of the ICR channels: 



A 
X 

w 



max ^ lo £2 (1 + 7fc) ( ! 8a ) 



fc=i 



s.t. Tr (NWW ff N fl ) < £ p (18b) 

Tr (WW") < P. (18c) 

The above problem is novel since the co-channel ICR interference term in the SINR is a 
function of the transmit signal itself. This is in sharp contrast with conventional single-cell 
E71 . multi-cell ESI , and underlay-only |fTT| | - [fT5ll downlink beamforming problems where the 
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co-channel interference is inevitably modeled as independent noise. While signal-dependent 
interference is a well-studied problem in radar signal processing [24), in our case this dependence 
manifests itself in a much more complicated and non-linear fashion involving exponential terms. 
We are faced with a non-convex objective function with multiple non-linear constraints, and at 
this point an analytical solution for W therefore appears to be intractable. 

To solve the sum-rate maximization problem numerically, a gradient projection (GP) algorithm 
can be used, which will converge to at least a locally-optimal stationary point. To compute the 
gradient of the UCR sum rate, we define the leakage term 

L k)j = J2 l^ w ;| 2 + h (W) + 4, (19) 



and compute Vw (R 

V Wfc (R, 



where 



^hfh.w^-lh.w,! 2 ^) 



(W 




-Zh k h k w k g^- 

(Ll,m) 



(20) 



dl k (W) _ 
d\v k 


V^ r> 2 9P Di 


= 


2 <^ w * ("*) /V A < (i-r) 

[hi ri (air 1 



(21) 



(22) 

and the differential on the RHS of (|2TT) is taken with respect to the average detection probability 
in (flOl) which is computable at the UCT 

At the k th iteration of the GP process, the updated precoding matrix W( fc ) in the direction 
of the gradient computed above will likely not satisfy the UCT transmit power and PR in- 
terference temperature constraints. The projection step of the GP algorithm therefore projects 
the iterate W^ back onto the feasible constraint set Q, + (P,£ P ), defined as £l+(P, £ p ) — {W | 
Tr (WW fl ) < P, Tr (NWW^N 1 ) < £ p }. Nominally, this would be achieved by determining a 
feasible W G &+(P, £ p ) that is closest to W^ in terms of Frobenius norm, i.e., by minimizing 
the squared distance d 2 (W , W<'>) = Tr ((W - W^) H (W - W^)) with appropriate 
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constraints: 

min d 2 (W ,W (fc) ) (23a) 

Wo 

s.t. Tr (W W^) < P (23b) 

Tr (NW„W^N ff ) < £ p . (23c) 

However, instead of numerically solving the above problem, a potentially suboptimal but much 
simpler approach is to scale W( fc ' such that both (I23bl) and (I23cl) are satisfied. This approach 
is partly motivated by the observation that the solution to (123 al) cannot satisfy both constraints 
with equality for a general channel N^I, and one of the constraints is guaranteed to be an 
inequality anyway. 

A summary of the GP approach for underlay prescient sum rate maximization is shown in 
Algorithm IIII-AL where the step sizes s k and a k are chosen using well-defined criteria such as 
Armijo's rule OH Sec. 2.3]. 



Algorithm III-A.l Prescient Gradient Projection Method 



Initialization: 

Set iteration index k — 0. 

Initialize W^ = [wf wf . . . wgj. 
Main Loop: 

1. Calculate the gradient V w o) (R s ). 

2. Choose an appropriate step size Sfe. Let W'( fc ) = W^ fc) + SfcV W (*o (R s ) 

3. Let W^ be the projection of W'^ onto fi + (P, £ p ), where 

Q + (P,Q 4 {W | Tr (WW H ) < P,Tr (NWW% H ) < £ p }. 

4. Choose appropriate step size a k . Let W( fc+1) = W (fe) + a k {W^ k) - wf } ) 

5. k = k + 1. If ||vec (W (fc) - W (fc_1) ) || < e, stop; else go to step 1. 



B. Algorithm Based on Convex Optimization 

While the iterative algorithm described above returns at least a locally optimal prescient beam- 
forming matrix, it is desirable to investigate designs based on simpler optimization procedures. 
In this section, we investigate a suboptimal approach that maximizes the partial underlay SINR 
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accounting for intra-UCR and PR interference, while making a best-effort attempt to limit the 
expected ICR interference by ensuring a minimum level of signal power leakage to them. We 
first define the partial UCR SINR B k as 

&= J^" Wfc L 2 , k = l,...,K u , (24) 

where the ICR interference term in the denominator of (fT6l) is omitted. Then we pose the problem 
of maximizing the minimum partial UCR SINR subject to a set of constraints {r)i} i=1 on the 
total UCT signal power received by the ICRs, as follows: 

max min B k (25a) 

w k 

s.t.Tr (WW") < P (25b) 

Tr (WW H F? F 4 ) >TH,i = l,...,K (25c) 

Tr (NWW ff N ff ) < £ p (25d) 

This can be posed as a convex optimization problem as follows. Let 3 k = w fc w", V&. Ap- 
plying a change of variable and relaxing the rank-1 constraints on 3 k , we have the reformulation 



max t (26a) 

{ J fc}fc=i 

s.t. t( J2 Tr (hf hi J,) + of J - Tr (hf h, J fc ) < (26b) 

Tr ((E^i <* Jfe ) F ^ F >rk,i = l,---,K (26c) 

J^^Tr^^P (26d) 

TrfNf^JfeJN^j <£ p (26e) 

t > (26f) 

J fc ^0, fc = l,...,/^. (26g) 

In this case, however, dropping the rank constraints on { Ji} i= " still does not lead to a semidefinite 
program (SDP), since the K u underlay SNR inequality constraints in (I26bl) are non-linear due 
to the fact that t is a variable. Therefore, a two-stage solution strategy is required where the 
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outer-loop performs a one-dimensional bisection search over t, while the inner loop solves (|26al ) 
for a given value of t, if feasible ll32l . 



C. Combined Downlink and Multicast Beamforming 

As an alternative suboptimal algorithm, we present an approach with a semi-analytical ex- 
pression for W, motivated by the simple observation that the detection probability of the energy 
detector in ([8]) increases monotonically with the received SNR at the SUs for a given false alarm 
rate Pfa,i- Consider the following two extreme cases for the choice of W: 

• Disregard ICRs, focus only on UCRs: If the UCT disregards the presence of the ICRs and 
focuses only on its intended receivers, a suitable choice for W is the RCI precoder W CT 
given by (flTT) . 

• Disregard UCRs, focus only on ICRs: At this extreme, the UCT ignores its downlink users 
and focuses only on improving the signal strength at the ICRs (particularly those that 
could produce the most interference). This is similar to a MIMO multicast (MC) downlink 
scenario, where priority is given to certain key users. A reasonable choice for the transmit 
precoder in this case would maximize the weighted average of the SNRs at the ICRs: 



K 

x 
w 



W MC = arg max £) ^Nja^Ti (F t WW H Ff ) , (27) 



■2 



where the weight P^Nja^ measures the interference impact of the ith ICR at the UCRs. 
The solution to (1271) is given by the dominant singular vectors of F^ E 9 F^ scaled by \[P, 
where F$ = £\ F.; and S 5 is a diagonal matrix with entries PiNja^, i = 1, • • • , K. 
Given that the prescient precoding objective is to balance these two competing goals, a sensible 
approach would be to choose W as some linear combination of the solutions: 

W, = aW C i + (1 - a) W M c < a < 1 , (28) 

where the optimal value of a E [0, 1] that maximizes (fTSl can be found by a simple line search. 

IV. MULTI- ANTENNA UNDERLAY RECEIVERS 

In this section we extend the prescient downlink precoding paradigm to the case of multi- 
antenna UCRs with multiple data streams transmitted to each. For simplicity, assume that each 
UCR is equipped with r u antennas, although the proposed algorithms hold for unequal array 
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sizes as long as the total number of receive antennas does not exceed t u . The extension to the 
case where the UCT serves t u spatial streams regardless of the total number of receive antennas 
can be made using the coordinated beamforming approach 11331 , for example. The received signal 
at UCR k is now 

K u K 

y k = H fc W fc s Uifc + ^ HfeWjS^- + ^ FiVk,i8i,i + n fc (29) 

&k i=i 

where H^ G C ruXt " is the main channel, W& G C tuXlk is the beamforming matrix applied to 

signal s U; fc G C^ xl for user k, F{ is the ICR indicator function as before, s/j is the i th ICR 

signal over interfering channel V k :i G C VuXri , and n k ~ CA/"(0, erf I) is additive Gaussian noise. 

The transmit covariance matrix for each UCR is given by Q& = W fc W^. We adopt a prescient 

block-diagonalization (PBD) strategy on the underlay downlink OBI . 041 to completely eliminate 

intra-UCR interference, as shown below. 

In the first approach, the transmit covariance matrices {Qk}f=i are computed jointly so as to 

optimize the underlay system sum rate while subject to constraints on the PR interference and 

the minimum power leaked to the ICRs. The proposed PBD scheme is described mathematically 

as 

max V log 2 1 1 + H fc Q fe Hf I (30a) 






s.t. HfcQjHf = 0,Vfc^j (30b) 



Tr ( N (Z) fc =! Q fc ) N ") ^ Zp (30c) 

Tr ( F > (EL"i Qfe ) F ^) " ^ Vz (30d) 

Tr (E^i <**) ^ P (30e) 

Q fc ^0, k = l,...,K u . (30f) 

Note that this is not equivalent to direct maximization of the UCR sum rate since the ICR 
interference is not included in the objective function. However, this decoupling renders the 
problem convex since the objective function is jointly concave and all constraints are linear in 
{Q fc }, and the leakage constraints r\i can be adjusted appropriately to diminish the probability 
of missed detections at the ICRs. 
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As an alternative PBD strategy, it is possible to separately design the precoding and power 
allocation matrices per user in a two-step process. Let 



H_fc 



Hi H fc _i H fc+1 H 



A'„ 



represent the overall UCR downlink channel excluding the k th user. First, a closed-form solution 
for the unit-power precoding matrix of user k is obtained from the nullspace of H_ fc . For example, 



from the SVD H_ fc = U_ fe S_ 



H 



the last (t u — Ik) right singular vectors 



contained in V_fc,o can be used to construct W& |33). However, unlike the conventional BD 
algorithm, the power allocated over the Ik spatial modes of user k is now no longer obtained via 
waterfilling. Let rank(HfcWfc) = r k for user fc's effective channel, and assume l k = r k - Consider 
the SVD of user k's effective channel H fc W fc = UfcE fc V^ where S fc = diag (e^i, . . . , efc,rj is a 
r k x r k diagonal matrix, and define A k = diag (X k ,i, . . . , Xk,r k ) to be the power allocation matrix. 
The overall downlink power allocation matrix is therefore A u = blkdiag (Ai, . . . , A^J. The PR 
interference and ICR signal power constraints are accommodated in the power allocation step 
based on a numerical optimization: 

K u l k 

™ x Yl Yl lo S2 (1 + 4,m x k, m ) (31a) 

S ' t - 5Z P 1 yiu "i 51 i H n c W fe,m|l2 X k,m > Cp (31b) 

*— »c=l *— 'k=1 *— 'm=l 

EV, V, ||fi,nW fc , m ||2 A fc , m > r/j,z = 1, ...,A", (31c) 

71=1 *— '/C=l ^— '771=1 

, > w fe>m LA fe < P (3 Id) 

fc=l * — 'm=l z 

where n c is the c th row of N, f in is the n" 1 row of Fj, and Wk, m is the m th column of Wfc. 
The leakage and power constraints (I31bl) - (l31dl) are equivalent to (I30cl) - (l30fl) . This is a convex 
program since the objective function is concave and all constraints are linear in {A^}, and can be 
efficiently solved using interior-point methods. It must be noted however that a separate design 
of the underlay precoding and power allocation matrices is potentially suboptimal compared to 
the joint design of (l30l) . 

V. Simulation Results 

In this section, we present the results of several numerical experiments to verify the im- 
provement in primary link performance provided by prescient beamforming. To avoid repetition, 
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unless specified otherwise, all results in this section are based on the partial CSI model with 
instantaneous CSI of the downlink and UCT-ICR links, and only statistical CSI of the ICR-to- 
UCR links available to the underlay transmitter. Each channel realization for all terminals is 
drawn from a zero-mean circularly symmetric complex Gaussian distribution, and all results are 
averaged over 1000 channel realizations. The background AWGN variance at all receivers is 
assumed to be unity, the primary antenna array sizes are fixed as t p = r p = 4 with PT transmit 
power P t = 10, and the PR interference cap is set to £ p = 10. The convex programs are solved 
numerically using the cvx MATLAB toolbox 11351 . At the ICRs we set the transmit power to 
Pi = 20dB, false alarm rate target Pfa,% = 10~ 3 Vi, and sample size of M = 4. The prescient GP 
algorithm is run 5 times for each set of channel realizations with four random initializations and 
an initialization based on the naive RCI precoder to reduce the likelihood of a local maximum; 
the best-performing precoding solution is chosen as the result. 

In Fig. 12, we first examine the energy detection receiver-operating-characteristic at an arbitrary 
ICR for prescient GP precoding compared to RCI transmission for K u = 3 single- antenna UCRs. 
The UCT transmit power is fixed at P = lOdB with t u — 3 antennas, and K — 2 ICRs are 
present with 77 = 2 antennas each. We observe that prescient precoding provides a significant 
improvement in energy detection performance for the entire range of Pfa, and consequently 
reduces the likelihood of ICR missed detections. 

Sum rate results for the single-antenna UCR downlink versus UCT transmit power with t u = 
K u = K = 3, Ti = 2 are shown in Fig. [3] The prescient schemes improve markedly upon 
the naive RCI precoder since each ICR with a missed detection interferes with multiple UCRs. 
The linear combination scheme is observed to be a very competitive alternative compared to 
the computationally intensive GP solution. The SDP-based prescient scheme suffers from the 
difficulty of optimally choosing leakage power thresholds 77$. The proposed prescient GP precoder 
provides an increase of up to 7 (bits/s/Hz) in spectral efficiency compared to the RCI scheme, 
which highlights the significant benefit of preemptively mitigating secondary user interference. 

We now consider prescient versus conventional block-diagonalization schemes for the multi- 
antenna UCR downlink with r u = 2. In Fig. [4] the underlay sum rate is displayed as a function 
of the UCT transmit power for t u — 8, K u = 4, K = 77 = 2. The greatest benefit of the PBD 
schemes is observed at low to intermediate SNRs, while the sum rate of all three algorithms 
gradually converge at high SNR. This is because the diversion of transmit power to the ICRs 
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under PBD now has a greater penalty in terms of spatial multiplexing loss to the multi-antenna 
UCRs. 

Finally, Fig. \5\ presents the PBD and BD underlay sum rates as the number of potentially 
interfering ICRs increases, for fixed UCT power P = 15dB. The relatively low combined 
transmit power of the UCT and PT leads to a potentially significant number of missed detections 
at the ICRs, and the expected ICR interference clearly worsens as K increases. This is especially 
true for the conventional BD scheme, which suffers from a pronounced degradation in sum rate 
since it neglects the sensing performance of the ICRs. An important implication of this outcome 
is that the successful coexistence of UCRs and ICRs in a heterogeneous DSA network cannot 
be assured merely by modifying the UCT precoding strategy; smarter ICR spectrum sensing 
approaches must also be adopted. 

VI. Conclusion 

This work has examined a novel heterogeneous DSA network where the primary users coexist 
with both underlay and interweave cognitive radios, all terminals being potentially equipped with 
multiple antennas. We investigated the design of MIMO precoding algorithms and the underlay 
transmitter in order to increase the detection probability at the ICRs, while simultaneously 
meeting a desired Quality-of-Service target for the underlay receivers and constraining the amount 
of interference leaked to the PUs. The objective of such a proactive approach, referred to as 
prescient precoding, is to minimize the probability of interference from ICRs to the UCR and 
PU receivers due to imperfect spectrum sensing. We presented three different downlink prescient 
precoding algorithms for the case of multiple single-antenna UCRs and multi-antenna PUs/ICRs. 
We then presented prescient block-diagonalization algorithms for the MIMO underlay downlink 
where spatial multiplexing is performed for multiple multi-antenna UCR receivers. Numerical 
experiments demonstrate that prescient precoding by the UCT preemptively mitigates missed 
detections at the ICRs, and provides a significantly pronounced performance gain in underlay 
sum rate compared to conventional precoding strategies. For future study, it is of interest to 
design prescient user selection algorithms to accommodate scenarios where the number of UCRs 
exceeds the number of UCT transmit antennas. 
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Fig. 1. Cognitive radio network with a multi-antenna underlay transmitter, K u underlay receivers, a single MIMO primary 
receiver, and K spectrum-sensing multi-antenna interweave cognitive radios. The primary transmitter and ICR-to-PR interfering 
links are not shown for clarity. 



21 



1 


1 1 1 — 


-B" 




- B - Prescient GP 


0.9 


- - * - RCI 




o 0.8 




"-i— « 


J2 " 


o 


* <" 


CD 


*■ y 


-1— « 

CD 

Q 0.7 




M— 


* ^ 


o 


' •* 


-1— « 


m ' 


1 0.6 


i y 


CD 
.O 


> ,' 


o 


i * 


» °- 5 


it 


CO 


i / 


a3 


a ' 


> 0.4 


i / 




• 


0.3f 


i 


i 


i i i i i i i i i 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Probability of False Alarm P p . 



Fig. 2. ROC curve for energy detection comparing prescient precoding with regularized channel inversion, t u = K u = 3, n 
K = 2,P= 1MB. 
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Fig. 3. Underlay sum rate for prescient algorithms and RCI precoding with t u = K u — K = 3, rj = 2. 
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Fig. 4. Underlay sum rate for prescient and conventional block-diagonalization with t u = 8, K u — 4, K — ri — r u — 2. 
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Fig. 5. Underlay sum rate for prescient and conventional block-diagonalization versus number of ICRs K, for t u = 6, K u 
3,r/ = r u =2,P= 100. 



